SEMIDEFINITE PROGRAMMING FOR PERMUTATION CODES 

MATHIEU BOGAERTS AND PETER DUKES 



Abstract. We initiate study of the Terwilliger algebra and related semidefi- 
nite programming techniques for the conjugacy scheme of the symmetric group 
Sym(n). In particular, we compute orbits of ordered pairs on Sym(n) acted upon 
by conjugation, explore a block diagonalization of the associated algebra, and 
obtain improved upper bounds on the size M(n, d) of permutation codes of 
lengths up to 7. For instance, these techniques detect the nonexistence of the 
projective plane of order six via M(6, 5) < 30 and yield a new best bound 
M(7i 4) < 535 for a challenging open case. Each of these represents an im- 
provement on earlier Delsarte linear programming results. 



1. Introduction and notation 

Consider the symmetric group Sym(n) on [n] := {1, . . . , n}. Let D be a collection 
of conjugacy classes of Sym(n). A subset T C Sym(n) is an (n, D)- permutation code 
if, for any two distinct elements (f>, tp £ T, their quotient cpip^ 1 belongs to a class 
in D. Note that the order of the terms is not important, since permutations are 
conjugate with their inverses. 

The Hamming distance between d) and ip in Sym(?i), denoted df/(0, is the 
number of non-fixed points of (j)ip . Equivalently, if and ip are written in single- 
line notation, df/(</>, "0) counts the number of disagreements between corresponding 
positions. It follows that dn is a metric for Sym(n). To the best of our knowledge, 
this metric was first considered by Farahat in [9]. 

In most investigations of permutation codes, the set D of admissible conjugacy 
classes reflects this metric and is instead taken as some subset D C [n]. This 
interpretation for D is that all conjugacy classes with n — d one-cycles, d 6 D, are 
allowed for quotients <fnp . In other words, in this context, an (n, £))-permutation 
code r is a subset of Sym(n) with all nonzero Hamming distances belonging to D. 
Based on the coding applications discussed below, the standard choice for D is the 
interval {d, d + 1, . . . , n} for some minimum distance d. 

Let M (n, D) denote the maximum size of an (n, Z))-permutation code. When 
D = {d, d + 1, . . . , n}, this is simply written M(n, d). Early work determining some 
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bounds and values for M(n, D) began in the late 1970s in Qj], where the term 
'permutation array' was used. Around 2000, permutation codes enjoyed a revival 
of interest with the discovery in [U [TU] of their applications to trellis codes. Then, 
the survey article [3] observed connections with permutation polynomials and also 
initiated the first serious computational attack on lower bounds on M(n,d). A 
probabilistic lower bound appears in |15j . Meanwhile, linear programming (LP) 
upper bounds were investigated, first in j!7j , and then subsequently in [21 [8] . A 
growing body of related work has emerged, including constant composition codes, 
injection codes, and study of the packing and covering radii. 

The LP bound rests on Delsarte's theory of association schemes [6] applied to 
the conjugacy scheme of Sym(n). More details follow in Section [5J In various 
combinatorial settings (block designs and binary codes, for instance) Delsarte LP 
bounds have been successfully improved using semidefinite programming (SDP). In 
general, the sizes of matrices required for this SDP tend to grow impractically large. 
So the most successful applications of SDP to designs and codes usually begin with 
an attack on the algebraic structure. In a little more detail, block diagonalizations 
of certain matrix algebras are desired in order to scale the computations. This is 
the content of Dion Gijswijt's dissertation [12]. 

For a finite set X, we use M X (C) or C XxX to denote the algebra of \X\ x \X\ 
matrices with rows and columns indexed by X . (Some canonical ordering of X is 
usually assumed, leading to the more standard notation M n (C) or C nx ", which 
we also use.) The (x,y)-entry of matrix A G M^(C) is here denoted A(x,y). All 
linear spans and generated sub-algebras (A±, A 2 , . . . ) are over C, unless otherwise 
noted. The conjugate transpose A T of A is denoted A* and as usual A is Hermitian 
if A* = A. A sub-C-algebra of Mx(C) closed under this Hermitian conjugate is 
called a C- (matrix) algebra. 



2. Bose-Mesner and Terwilliger algebras 

An d- class association scheme on a finite set A is a list of binary relations 
Ro, . . . , Rd on X such that Ro is the identity relation, the relations partition X 2 , 
and the following regularity condition holds: given x and y with (x, y) G Rk, 
the number of z G X for which both (x,z) G Ri and (z,y) G Rj is a constant 
depending only on i,j,k G {0, ...,d}. These values p^ are called the structure 
constants or intersection numbers of the scheme. Here, we are mainly interested 
in symmetric schemes in which each relation Ri is symmetric and thus p\j = p^. 
Chris Godsil's notes [13] offer a readily available and comprehensive reference on 
association schemes. 

Example 2.1. Consider X — G, a finite group. The relations Ri are indexed by 
the conjugacy classes of G, where i? corresponds with the trivial conjugacy class. 
Declare (g, h) in Ri if and only if ghr 1 belongs to the ith conjugacy class. The fact 
that these Ri partition G 2 is clear. A character sum can compute the structure 
constants (see below). This is the conjugacy scheme on G. 
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Let \X\ = N. Define the N x N zero-one matrices Ao, . . . , Ad G Mx(C) by 
A(x,y) 



I if (x, y) G i?i 
otherwise. 



We have Ao = J and Ao + • • • + Ad = J, the all ones matrix. And from the definition 
of the structure constants, 



(2-1) AA^Y.Pij^ 



k=0 

From this it follows that (Aq, . . . , Ad) is a commutative C*-matrix sub-algebra of 
Mx(C) of dimension d + 1. This is the Bose-Mesner algebra 21 of the scheme. 

From spectral theory, 21 also has a basis of orthogonal idempotents Eq , . . . , E^ 
with Eq + ■ ■ ■ + Ed = I ■ It is convenient to index Eq = j^J, which is evidently one 
of the idempotents. 

The Hadamard or entrywise product of matrices A, B G <C NxN is A o _B with 
(A o B)(i,j) = A(i,j)B(i,j). For the Bose-Mesner algebra, we have the dual 
relations 

Ai o A.j = SijAi and /•.', /•/, = SijE,. 
In particular, since 21 is closed under o, a parallel version of (12 . If) exists; namely, 

1 d 

fe=0 

for some positive reals . These are called the Krein parameters for the scheme. 

The basis change matrices between (Ao, . . . ,Ad) and (Eo, ■ ■ ■ ,Ed) are denoted 
by P and jjQ', explicitly, the entries are given by 

d 1 ^ 

Ai = J2 P(hj)Ej and Ei = - £ Q(i, j)A,, 

3=0 j=0 

The -Ej being projections imply that the entries P(i,j) are the eigenvalues of the 
matrices Aj, with the column space of Ej as the associated eigenspace. For this 
reason, the matrices P, Q are called, respectively, the first and second eigenmatrix. 
They are of course related by Q — NP^ 1 . 

For a subset 7 CI, its indicator vector is ly G {0, 1} , with 1y{%) = 1 if and 
only if x G Y . The inner distribution of a nonempty Y is a = (ao, . . . , ad), where 



The following observation leads to the famous Delsarte LP bound for cliques in 
association schemes. 

Theorem 2.2 ((Delsarte, [B])). The inner distribution & of a nonempty subset of 
points in an association scheme with second eigenmatrix Q satisfies aQ > 0. 
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Now consider a pointed set (X, x) with x <E X. Define diagonal zero-one matrices 
E' ,...,E' d mM x (C)by 



1 if (x, y) 6 Ri 
otherwise. 



Observe that, by analogy with the Bose-Mesner idempotents, these are projections 
with Eq of rank one and E' + ■ ■ ■ + E' d = I. 

The algebra 1 = (AijE'j) obtained by extending 21 by these E[ is called the 
subconstituent or Terwilliger algebra of the scheme with respect to x. Unlike 21, the 
algebra T is not in general commutative. 

Completing the duality, we have A' Ql . . . , A' d with 

A' i {y,y) = NE i {x,y). 

After this introduction, we are now interested exclusively in the conjugacy scheme 
on Sym(n). Recall that the conjugacy classes (and irreducible characters) of Sym(ri) 
are in correspondence with the integer partitions of n. For a fixed n, we will use m to 
denote the number of partitions, and index the conjugacy classes as Co, . . . , C m _i, 
with Co = {id} consisting only of the identity permutation. For any fixed 6 E Cfc, 
we have 

tf !i = \{(<t>,il>)eC i xC j :<fnl> = 0}\. 
Alternatively, as stated in [14] . 

k _ ICHQI v x(^)x(0 3 )x(^) 
PlJ n\ ^ X (id) 

where the sum is over all m irreducible characters, and where fa is a representative 
in class i, etc. 

The Bose-Mesner algebra 21 is a commutative sub-algebra of M n \(C). Its eigen- 
values are obtained nearly directly from the character table of Sym(n). 

Lemma 2.3 (Tarnanen, [17]). The second eigenvalue matrix Q for the conjugacy 
scheme on Sym(n) is given by 

Q{hi) = Xj(id)Xj(<l>i), 
where < i,j < m index both the conjgacy classes and irreducible characters. 

The idempotents and Krein parameters of this conjugacy scheme are also ex- 
pressible via the irreducible characters, but we omit the details. 



The linear programming bound of Theorem 12.21 now becomes 

maximize: oq + a\ + ■ ■ ■ + a m _i 

subject to: J2o<i <m a iXk((f>i) > for < k < m, 

ao = 1, a,i > 0, and 

a, = if dff(id, (pi) <j£ D. 
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Given n and D, let Mlp denote this maximum value. Then 

M(n, D) < M LP . 

This was used with some success in [31 [T7] to obtain some upper bounds on permu- 
tation codes with n < 15. Additionally, the LP led to a general upper bound on 
M(n,4) in 0. 

When extending to the Terwilliger algebra in this setting, it is natural to take 
id G Sym(n) as the distinguised element. We have 



l if e Ci, 

otherwise. 



The products E'fiLE'j are supported on the Cj x Cj block. This offers a natural 
decomposition of the Terwilliger algebra in general. Note that E^A^E'^ is the zero 
matrix if and only if = 0. 



3. ISOMETRIES, ORBITS, AND THE CENTRALIZER ALGEBRA 



Consider the automorphism group Aut(2t) of the Bose-Mesner algebra of Sym(n). 
On one hand, this is the intersection of the automorphism groups of the graphs 
whose adjacency matrices are the zero-one generators A4, < i < m = p(n). 
Alternatively, this is also the group of isometries of Sym(n) , when endowed with the 
Hamming distance metric dn- Let us also denote this group by Iso(n). 

Let C and 1Z denote the group actions of left and right multiplication by Sym(n), 
acting on itself. Let I be the (involutoric) group action generated by inversion 
<f) h-> (j)^ 1 on Sym(n). That these actions induce all isometries appears in [5], and is 
discussed further in [2]. 

Lemma 3.1. Aut(2t) = Iso(n) = (£ x K) x 1 = Sym(n) I 2. 



Let Isoi(n) be the subgroup of isometrics which fixes the identity element. Its 
action on Sym(n) is generated by conjugations C and inversion X. These actions 
commute, and therefore Isoi(n) = C xT = Sym(n) x 2. 

For the extension to the Terwilliger algebra T, we are concerned with how Isoi (n) 
acts on pairs of permutations. In what follows, sums indexed by 'x' are to be taken 
over all irreducible characters \ °f Sym(n). 

Proposition 3.2. The number of orbits o/Sym(n) 2 , acted on by conjugation and 
inversion is given by 



aGSym(n) 
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Proof. We use Burnside's orbit counting lemma to obtain 

6 « = ]cVT\ £ |fix(5)l 



gGCxX 



1 

2d 

1 

2d 



53 :a(0,V)a- 1 = (^»}|+ 53 |{(& V) : o(^, ^)o _1 = (<T \ lT X )}| 

aGSym(n) aESym(n) 



53 \{cf> : a^a' 1 = 0}| 2 + ^ K^ra^a- 1 = 

aGSym(n) aESym(n) 



Both terms in the sum are connected with interesting class functions on Sym(n). 
The first term is just the sum of squared centralizers. Recall that orthonormality 
of the irreducible characters implies that for an element a 6 Sym(n), 



53x(a) 2 = |C(a)|, 



the (size of the) centralizer of a. So the first of our two summation terms can be 
rewritten 



(3.1) 



QGSym(n) \ X / 



For the second sum, observe that the condition a.§a~ x = _1 is equivalent to 
(a(j)) 2 = a 2 . So for a fixed a, the number of such <j> is simply the number of square 
roots of a 2 in Sym(n). As an aside, this counting problem for general groups G was 
connected long ago (see [18]) with the 'Frobenius-Schur indicator' 

s (x) = ]77| £ 
11 sec 

For our purposes, all irreducible representations of Sym(rt) are real and therefore 
s(x) is always 1. The number of square roots of an element j3 reduces in this case 
to the basic character sum ^ xifi)- For j3 = a 2 , the term we seek becomes 



E £*(< 

aGSym(ri) \ X 



= E (ex^VeW) 

,36Sym(ri) \ X / \ X / 



(3.2) 



= E £*(< 

QGSym(n) \ X 



Combining (13. ip and (|3.2[) completes the proof. 



□ 



Now consider IsoJ(n) = Sym(n) x 2 x 2, acting on Sym(n) 2 by conjugation, 
inversion, and the 'coordinate swaps' (4>,ip) l— * (V'i^)- A small extension of Propo- 
sition 13.21 also yields a character sum formula to count orbits for this action. 
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Proposition 3.3. The number of orbits ofSym(n) 2 , acted on by conjugation, in- 
version, and coordinate swaps, is given by 



where b n is as in Provosition \3.2\ 



aSSym(ra) V X / \ X / 



Proof. New sums not already present in the proof of Proposition 3.2 count, for a 
fixed a, 

|{(&V0 : a(0,V)ar x = (ip,4>)}\ and |{(0, V) : afo^e*" 1 = (i/r 1 ,^- 1 )}!- 

Since conjugacy by a is a bijection, it follows that each of these is counted simply by 
the centralizer of a 2 . Alternatively, we may count, over all (3 6 Sym(n), the number 
of square roots of j3 times the size of the centralizer of j3. Using the expressions in 
the previous proof, 

E i<V)i= E fexw) fe*W). 

aeSym(n) /3GSym(n) V X / V X / 

This, along with earlier terms, finishes the orbit count. □ 



Table [T] contains some small values of b n and 6* . 

Table 1. Orbits of Sym(n) 2 under conjugation with inversion, and 
with inversion and swaps. 



n 


n\ 


b n 


b* 


4 


24 


43 


28 


5 


120 


155 


93 


6 


720 


761 


425 


7 


5040 


4043 


2151 


8 


40320 


27190 


14016 



Let Oi, . . . , Om denote the orbits of Sym(n) 2 under Isoi(n) = Sym(n) x 2. Define 
the corresponding zero-one matrices B\, . , . , Bm G M n \(C) by 

I otherwise. 

Put 03 = (Bi : I = 1, . . . , M) and observe that by Proposition l3.2l we have dim(03) = 

b n - 

Alternatively, 03 = Endc(A), where X = Sym(n) and G — Isoi(n). For a 6 
Sym(n), consider the zero-one matrix a G M n \(C) with 

_ jl if a~ x ^a = if), 

I otherwise. 
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Then a *-y a is the conjugacy representation. Define the action of inversion similarly 



Let A = (i, a : a £ Sym(n)). The following rephrases the earlier definition of orbits 
under Isoi(n) in terms of matrix algebras. 

Proposition 3.4. 03 is the centralizer algebra of R. 

Since conjugacy classes of Sym(n) are left invariant by both conjugation and 
inversion, the orbits under Isoi(n) admit a classification according to conjugacy 
classes. In other words, each matrix Bi is supported on some block x Cj. The 
sum of all such Bi supported on a given block is, when restricted to that block, the 
all-ones matrix J. 

With notation from Section 2, we recall that the 'restricted' Bose-Mesner algebra 
E'jI&Ej has conjugation-invariant generators supported on the d x Cj block. So it 
follows that T is a sub-C-algebra of 03. The reverse inclusion seems reasonable, but 
much more difficult to prove. Indeed, the dimensions in Table [T] grow much faster 
than the number of generators for T. 

Conjecture 3.5. For the conjugacy scheme o/Sym(n), the Terwilliger algebra has 
as a basis the zero-one matrices Bi, i = 1,...,M, defined via orbits of Sym(n) 2 
under Isoi(n). That is, as C-algebras, we have T = 03. 

Conjecture l3.5l was verified up to n — 5 in pQ . We have extended this verification 
to n — 6 but leave the general problem to future work. 

Fortunately, the center of 03 is reasonably easy to describle. The following 
appears in |16j for general groups; here (and from now on) we are interested in 
G = Isoi(n) = Sym(n) x 2. 

Proposition 3.6. The primitive central idempotents of 03 are given by 



where < k < 2m indexes the irreducible characters of G. 

Since the commute, they have a common basis of eigenvectors. Consider a 
unitary matrix U G M n \ (C) which has orthonormal eigenvectors of the e^, organized 
by columns. Then U ( 3U T is block-diagonalized by simple blocks. Our goal later, 
in Section 5, is to decompose these simple blocks into basic blocks. For now, 
though, we have sufficient background to introduce semidefinite programming for 
permutation codes. 




(3.3) 




(id) 



7£G 



4. Semidefinite programming 



Let r be a subset of Sym(n), and recall that lr is the n! x I zero-one indicator 
vector of T. 
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Define now the subgroup II of Iso(n) consisting of all isometries t such that id £ 
i(r). Let II' denote the complement of this subgroup, consisting of all t' € Iso(n) 
with id £ t'(T). It is clear that |II| = 2n!|r|, while |n'| = 2n!(n! - |r|). This defines 
two matrices 

Rt = lff| X] 1 *( r ) ■ 1 *(r) and R 'r = ijtTT $Z 1 *'(r) ■ ^(r) 
1 1 ten 1 1 ferr 

which are roughly analogous to the inner distribution vector a of Section 2. 

By construction, both Rr and R T are positive semidefinite (symmetric) matrices. 
And the trace of Rr gives the cardinality of T. The (id, id)-entry of Rr is always 1, 
since id appears in all copies under II. The the (id, 0)-entry equals the ((/>, </))-entry 
in Rr, and, Rr being symmetric, it is also equal to the (</), id)-entry. Since Iso(n) 
is transitive on Sym(n), one has the following formula and lemma. 

(4.1) flr(6l/0 + ! ^^rW>,V>) = flrf^f^f 1 ) = iZrO'W 1 )- 

(In particular, the matrix i?r + " ! | rP -R r has a diagonal of 1.) 

Lemma 4.1. Consider an entry (<f>,ip) S Sym(n) 2 . //no quotient in T belongs to 
the same conjugacy class as cftip' 1 , then i?r(</>, VO = R r (<j),ip) = 0. 

Next, we have the key connection with the centralizer algebra 03 discussed in 
Section 3. 

Proposition 4.2. The two matrices Rr and R' T belong to the algebra 03. 



Proof. Recall that {/?/} are the zero-one matrices representing the orbits {O/} of 
Sym(n) 2 under the action of Isoi(n). With the convention that O-y = {(id, id)}, 
one has Rr = B± + a/_B/, where the coefficients a; are explicitly determined by 
the formula 

_ \{(e,fi,v) g r x r x r : (fl-y,^- 1 ^) e o,}\ 

\T\\Oi\ 

Observe that < ai < 1. 

Let 1\ denote the subset of indices I ^ 1 such that Oi contains elements of type 
(id, </>), (cf), id) or (</>, </>), and let I2 be the complementary set {2 . . . M}\I\. 

With some rearranging and (|4.1I) . it follows easily that 
n!-|r| 



l -R' r = (I-B 1 )+J2aiB' l , 



1 1 1 
where the matrices B[ are defined by B[ = —Bi if I £ I2 and by the following 
expression if Z £ I\ : 

( 1 if (^-\id) eO h 

B'li^ i>) = \ -1 if <t> = i> an d (<f>, id) £ O u 
y otherwise. 

In both cases, we have the required matrices expressed as linear combinations of 
the {Bi}. □ 



10 



MATHIEU BOGAERTS AND PETER DUKES 



Wc now summarizic the above facts concerning Rr and R' r . 

Proposition 4.3. Let D C [0,n]. For any (n,D) -permutation code T C Sym(n), 
we have 

' |r|-Tr(i? r ) 
Rr = B 1 + £). ai Bi h 

R^=I-B 1 +'£.a j B f j hO 
< # r (id,id) = 1 
i?^(id, id) = 

R T {<j>^) = R' T { ( j>,<il>)=Q ifd H (<f>,il>)tD 
{ i?D 

From this, the following SDP system is obtained, considering the coefficients cij 
as variables in [0, 1]: 

maximize: Tr(i?i) 

subject to: i?i = B\ + ^ aiBi >; 

Ri = I - B 1 + £\ aiB'i h0, 0<a t <l 

Hi (id, id) = 1, 

fl 2 (id, id) = 0, 

R 1 (<P,ij) = R 2 (cj>,^)=0 if dui^) <£D. 

Given n and D, let Ms dp denote the maximum value of this program. Then 

M(n,D) < Msdp- 



Table 2. LP and SDP bounds for Sym(6) 



D 


M S dp 


M LP 




D 


M S dp 


M LP 







1 


1 


1 


{4,5,6} 


120 


120 


120 


{6} 


6 


6 


6 


{3,5,6} 


56 


60 


90 


{5} 


12 


15 


5 


{2,5,6} 


27 


30 


60 


{4} 


7 


12 


4 


{3,4,6} 


39 


48 


72 


{3} 


5 


6 


3 


{2,4,6} 


48 


48 


48 


{2} 


2 


2 


2 


{2,3,6} 


20 


36 


36 


{5,6} 


25 


30 


30 


{3,4,5} 


60 


60 


60 


{4,6} 


35 


48 


24 


{2,4,5} 


24 


30 


80 


{3,6} 


18 


24 


18 


{2,3,5} 


14 


15 


30 


{2,6} 


9 


12 


12 


{2,3,4} 


24 


24 


24 


{4,5} 


20 


20 


20 


{3,4,5,6} 


360 


360 


360 


{3,5} 


14 


15 


15 


{2,4,5,6} 


120 


120 


240 


{2,5} 


12 


15 


10 


{2,3,5,6} 


56 


60 


180 


{3,4} 


15 


24 


12 


{2,3,4,6} 


48 


48 


144 


{2,4} 


8 


12 


8 


{2,3,4,5} 


120 


120 


120 


{2,3} 


6 


6 


6 


{2,3,4,5,6} 


720 


720 


720 



Recall that the matrices Bi have size n! x n\. So presently, we can only directly 
consider this SDP problem for n < 6. Since permutation codes are well understood 
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for n < 5, our main contribution at this stage is a table of SDP bounds for n = 
6. For n > 7, it is necessary to consider an equivalent SDP obtained via block- 
diagonalization of the algebra 03. That is the topic of the next section. 

We conclude with Table 2, indicating the SDP bounds for n = 6 and various 
distance sets D. For comparison, the Delsarte LP bounds M^p are given, along 
with the quantity lideDd, which is an upper bound cither for M(n, D) or M(n, D c ). 



First, let us recall the notion of ^-isomorphism. A linear mapping n : A i— > B 
such that for all (A, B) 6 A x B, k(AB) = k{A)k(B) and k(A*) = n(A)* is called a 
*-homomorphism of algebras. A bijective such mapping is called a C* -isomorphism 
and the two algebras A, B are *-isomophic. A ^-isomorphism preserves positive 
semidefiniteness. 

We may regard & as a semisimple G-module, acted upon by *B. It therefore de- 
composes as a sum of orthogonal submodules Wk , where Wk has columns forming 
a basis of the eigenspace of the matrix e^, defined in (13. 31) . for < k < 2m. Each 
Wk is a direct sum of m k submodules: Wk = V k x © ■ ■ ■ © V k " lk , with V k l ~ V£ for all 
Let dk denote the dimension of Vi- 

Now we would like to decompose each submodule W k of 03 = Endc(^) = © fc Wk 
into smaller blocks (that we call basic blocks). The method is the following: for 
each k, compute an element Vk G Wk- Then define Uk '■= 23i>fc, that is a linear 
subspace of Wk ■ 

Note that the vector Vk can be obtained as ekX for any vector x, because 
ek{ekx) = CkX since e& is idempotent. So any non-zero row of tk would give an 
appropriate element of Uk- 

By Schur's lemma, the subspace Uk meets each VJ* in a 1-dimensional subspace. 
Then for all Bj € 03, Bj : Uk — > U k has an action on the basis of Uk- If ("i, ■ ■ ■ , u mk ) 
is an orthogonal basis of Uk , then the basic block Bik := Bi\u k of Bi can be obtained 



Each basic block Bik is an mk X matrix. The matrix Bi has nonzero entries in 
(<fi, ij}) position if ((fiyip) 6 Oi. The scalar product defining Bik can be written as: 



5. Block diagonalization 



by: 



(B ik )(j,l) 



\\ u j\\ \\ u i\\ 
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Table 3. LP and SDP bounds for Sym(7) 



D 


M S D p 


M LP 




D 


M s DP 


M LP 







1 


1 


1 


{4,5} 


26 


60 


20 


{7} 


7 


7 


7 


{3,5} 


15 


15 


15 


{6} 


28 


30 


6 


{2,5} 


14 


15 


15 


{5} 


14 


15 


5 


{3,4} 


21 


26 


12 


{4} 


9 


12 


4 


{2,4} 


9 


12 


8 


{3} 


6 


8 


3 


{2,3} 


7 


9 


6 


{2} 


2 


2 


2 


{5,6,7} 


134 


140 


210 


{6,7} 


42 


42 


42 


{4,6,7} 


181 


205 


168 


{5,7} 


47 


52 


35 


{3,6,7} 


70 


84 


126 


{4,7} 


34 


46 


28 


{2,6,7} 


59 


84 


84 


{3,7} 


27 


42 


21 


{4,5,7} 


74 


93 


140 


{2,7} 


12 


14 


14 


{3,5,7} 


54 


63 


105 


{5,6} 


30 


30 


30 


{2,5,7} 


47 


52 


70 


{4,6} 


47 


72 


24 


{3,4,7} 


49 


84 


84 


{3,6} 


31 


36 


18 


{2,4,7} 


36 


46 


56 


{2,6} 


34 


48 


12 


{2,3,7} 


32 


42 


42 



D 


M SDP 


M LP 


H-deod 


D 


M SD p 


M LP 


RdeDd 


{4,5,6} 


120 


120 


120 


{3,4,5,7} 


81 


93 


420 


{3,5,6} 


66 


108 


90 


{2,4,5,7} 


90 


140 


280 


{2,5,6} 


42 


60 


60 


{2,3,5,7} 


54 


63 


210 


{3,4,6} 


50 


72 


72 


{2,3,4,7} 


77 


168 


168 


{2,4,6} 


48 


72 


48 


{3,4,5,6} 


360 


360 


360 


{2,3,6} 


35 


54 


36 


{2,4,5,6} 


120 


120 


240 


{3,4,5} 


60 


60 


60 


{2,3,5,6} 


77 


108 


180 


{2,4,5} 


32 


60 


40 


{2,3,4,6} 


51 


72 


144 


{2,3,5} 


15 


15 


30 


{2,3,4,5} 


120 


120 


120 


{2,3,4} 


24 


32 


24 


{3,4,5,6,7} 


2520 


2520 


2520 


{4,5,6,7} 


535 


543 


840 


{2,4,5,6,7} 


612 


630 


1680 


{3,5,6,7} 


380 


420 


630 


{2,3,5,6,7} 


380 


420 


1260 


{2,5,6,7} 


160 


172 


420 


{2,3,4,6,7} 


207 


280 


1008 


{3,4,6,7} 


207 


280 


504 


{2,3,4,5,7} 


120 


168 


840 


{2,4,6,7} 


181 


205 


336 


{2,3,4,5,6} 


720 


720 


720 


{2,3,6,7} 


74 


84 


252 


{2,3,4,5,6,7} 


5040 


5040 


5040 



The algebra 03 = ©^q 1 C mfcXmfc generated by the diagonal joins of Bit yields 
an algebra that is *-isomorphic to 03. Let us dchne the linear mapping 







$ : Bi Di 





Bi2 



B, 



By construction, $ is a ^-isomorphism, inducing the equivalence between the two 
systems below. 



SDP FOR PERMUTATION CODES 
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maximize: Tr(X) ( maximize: Tr(X) 

subject to: X = ]T\ a l C l ^ <=>■ I subject to: $(X) = £V o»^(C») h 
< a 4 < 1 ( < a, < 1 

After similar computations as in the previous section, we can report SDP bounds 
for n = 7 according to Table 3. 



6. Concluding remarks 



For n > 8, the number of variables is presently too large for practical imple- 
mentation of the SDP, even after block-diagonalization. But we are hopeful that 
n = 8 can soon be attacked with some new ideas for generating orbits, block- 
diagonalizing, and possibly constraining the number of variables. In fact, should 
the SDP eventually scale up to n = 10, there is a chance that M(10, 9) < 90 can 
be proved in this way. Following [4 , this represents an opportunity for a possible 
alternate proof of the nonexistence of projective planes of order 10. 
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